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Abstract 

The properties of the ground state of 178 Hf and the isomeric state ^Hf 7 ™ 2 are studied within 
the adiabatic and diabatic constrained relativistic mean field (RMF) approaches. The RMF 
calculations reproduce well the binding energy and the deformation for the ground state of 
178 Hf. Using the ground state single-particle eigenvalues obtained in the present calculation, 
the lowest excitation configuration with K n = 16 + is found to be ^(7/2~[514]) _1 (9/2 + [624]) 1 
7r(7/2 + [404]) _1 (9/2 _ [514]) 1 . Its excitation energy calculated by the RMF theory with time-odd 
fields taken into account is equal to 2.801 MeV, i.e., close to the i'SHf"* 2 experimental excitation 
energy 2.446 MeV. The self-consistent procedure accounting for the time-odd component of the 
meson fields is the most important aspect of the present calculation. 
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The relativistic mean-field (RMF) theory is one of the most successful microscopic models 
in nuclear physics. &2>& From the very beginning, it incorporates the important relativistic 
effects, and it has achieved success in describing many nuclear phenomena related to stable 
nuclei exotic nuclei as well as supernova and neutron stars ®. The RMF theory 
provides a new explanation for the identical bands in superdeformed nuclei^ and for the 
neutron halo in heavy nuclei ^ , it predicts giant neutron halos, a new phenomenon in heavy 
nuclei close to the neutron drip line ^'■ 8 -', it naturally generates the spin-orbit potential, 
explains the origin of the pseudospin symmetry as a relativistic symmetry [E>i2»^i) ) and spin 
symmetry in the anti-nucleon spectrum and also describes well the magnetic rotation 
the collective multipole excitations ^ as well as the properties of hypernuclei^, etc. Lately, 
the ground state properties of about 7000 nuclei have been calculated in the RMF+BCS 
model and good agreements with existing experimental data were obtained Recent and 
more complete reviews of the applic ations of the RMF model, particularly, those to exotic 



nuclei, can be found in Refs. 
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Recently the 31-yr isomer of 178 Hf (also called ^Hf™ 2 , K* = 16+, E x =2AA6MeV) 
has attracted extensive attention ^ — '— ' — ' for its potential to be a good medium 
of energy storage «£L The long half-life of ^Hf™ 2 is connected with the strong in- 
hibition of spontaneous electromagnetic transitions restricted by the K selection rule, 
and it supports the point of view that this high-.fr state has an axially symmetric 
intrinsic shape The configuration originally suggested for the isomer 178 Hf n2 is 

z/ 2 (7/2-[514])(9/2 + [624])tt 2 (7/2 + [404])(9/2-[514]) f25 ' , which was further supported by the 
view of alignment and by the g-factor of the corresponding rotation bands. ^ 

In this Letter, we will study the properties of the ground state of 178 Hf and also investi- 
gate the possible configuration of ^Hf™ 2 within the self-consistent axially symmetric RMF 
theory. From the point of views of the adiabatic constrained calculation, nuclear ground 
state presents the global minimum of the potential energy surface (PES). To study the 
excited states, like isomers, the diabatic (configuration-fixed) constrained approach can be 
applied as an effective method. »ZJ Within the adiabatic constrained approach, nucleons al- 
ways occupy the lowest levels, while in the diabatic constrained approach, the configuration 
is kept fixed by the so-called concept of "parallel transport" . In the present paper, both the 
adiabatic and diabatic constrained RMF approaches are applied in our investigation. 

The basic ansatz of the RMF theory is a Lagrangian density where nucleons are described 
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as Dirac particles which interact via the exchange of various mesons and the photon. The 
mesons considered are the isoscalar-scalar a, the isoscalar-vector u and the isovector-vector 
p. The effective Lagrangian density reads «J 



C = ip 
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in which the field tensors for the vector mesons and the photon are, respectively, defined as 

(2) 

From the Lagrangian, the equation of motion for the nucleon is 

{a . [_ zV - V(r)} + V (r) + (3[M + S(r)]}^ = e^, (3) 

with the attractive scalar potential S(r) = g u u(r), the usual repulsive vector poten- 

[r) + e 
l-r 3 



1 - r 3 

tial V (r) = gu>u (r) + g p r 3 p (r) + e — - — A (r) , and the nuclear magnetic potential 



V r) 



9u,u[r) + g P T 3 p{r) 



-A(r). The Klein-Gordon equations for the mesons 



and electromagnetic fields are 



-V 2 + m^)C(r) = S c (r), 



(4) 



where S^(r) is the source term and all other notations are the same as in Ref. 
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In the RMF approaches which are widely used, only the time-even fields are essential for 
the physical observables, since the time-odd components of vector fields do not exist because 
of the time reversal symmetry for the ground state of an even-even nucleus. For an odd- A 
or odd-odd nucleus, the unpaired valence nucleon will give non-vanishing contribution to 
the nuclear current which provides the time-odd component of vector fields, i.e., the nuclear 
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magnetic potential. It is found that the nuclear magnetic potential has small influence on 
the root-mean-square radii and quadrupole moments while it plays an important role in 
the single-particle properties and magnetic moments in odd- A or odd-odd nuclei. One 
should keep in mind that for the excited states in the even-even nuclei, there may also exist 
the unpaired nucleons, which result in non- vanishing time-odd field. Therefore the time-odd 
fields should also be treated carefully for some isomeric states in the even-even nuclei, as 
in the odd- A or odd-odd nuclei. In the calculation without current, the nuclear magnetic 
potential V(r) will be neglected. 

For the adiabatic constrained approach, the binding energy at a certain deformation is 
obtained by constraining the mass quadruple moment (Q 2 ) to a given value //, i.e. 

(H') = (H)+ l -C{{Q 2 )-tf. (5) 

where C is the curvature constant parameter, and li is the given quadrupole moment. The 

expectation value of Q 2 is (Q 2 ) = (Q2)n + (Q2) P , where {Q 2 ) n ,p = (2r 2 P 2 (cos 9)) n , P - The 

3 

deformation parameter (3 2 is related to (Q 2 ) by (Q 2 ) = Ar 2 (3 2 , r = RqA 1 ^ 3 (R = 1.2 

V57T 

fm) and A is the mass number. By varying fi, the binding energy at different deformations 
can be obtained 

For the adiabatic constrained approach, the occupied levels are determined by the so- 
called "parallel transport" ^2, i.e., 

(^(9)1^(9 + Ag))| A9 -*o ~ Sij, (6) 

where i and j enumerate all the single-particle levels of two adjacent configurations. In 
such a way, the original configuration at q can be traced and the corresponding PES can 
be obtained as a function of the deformation ^- 7 h In principle, if Aq is small enough, the 
configurations at q and at q + Aq should be the same. In the calculation, the two-step 
procedure is adopted: first, the wave functions and the configuration at the initial q are 
recorded. Second, the wave functions \ipi(q)) are mapped to \ipj(q + Aq)) one by one by 
searching the largest overlap in \ipj(q + Aq)) with the same quantum number Q 71 . The 
configuration is transferred by copying the occupation number from \ipi{q)) to the mapped 
\ipj{q + Aq)). The wave functions and the configuration at this q + Aq are also recorded. 
The second step is repeated until enough points on the diabatic PES are obtained. 

The constrained RMF calculations are carried out with parameter set PK1 The full 
N = 20 deformed harmonic-oscillator shells for fermions and bosons are taken into account 
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as the basis. This basis is large enough to produce a converged binding energy at certain 
deformation. 

The PES of 178 Hf obtained in adiabatic (open circles) and diabatic (lines) constrained 
RMF calculations are plotted in Fig. [TJ The calculated energy £'=-1434.0 MeV and the 
deformation /3 2 =0.283 of the ground state (denoted as a black asterisk in Fig. [TJ are in good 
agreement with the experimental energy -1432.8 MeV t22Jand deformation 0.280 f33 L In this 
figure, the adiabatic PES can be decomposed into three regions by the discontinuity, i.e., 
/? 2 =0.22 ~ 0.32 (region 1), (3 2 = 0.35 ~ 0.40 (region 2), f3 = 0.43 ~ 0.50 (region 3). It is 
known that the discontinuity originates from the change of the configurations, and the config- 
uration for all points in one region is the same. This is confirmed by the diabatic constrained 
calculation in which the configuration is kept fixed during the constraint procedure. It can 
be seen in Fig. [1] that for every region the diabatic calculation (solid curves) coincides with 
the adiabatic calculation (open circles) and extends the region much wider. The crossover 
of the solid curves announces the change of the configurations. For example, from region 1 
to region 2, the configuration changes from the ground state to 7r(7/2 + [404])~ 2 (l/2~[541]) 2 , 
where microscopically the proton level 7/2 + [404] below the Fermi surface in region 1 be- 
comes unoccupied while the proton level 1/2~[541] above the Fermi surface in region 1 
becomes occupied. Since a pair of protons change the levels at the same time, the K-values 
remain zero for region 2. Similarly, from region 2 to region 3, the configuration changes 
from 7r(7/2+[404])- 2 (l/2-[541]) 2 to z/(5/2-[512])- 2 (1/2 + [660]) 2 tt(7/2 + [404])- 2 (1/2-[541]) 2 
and the K- values also remain zero for region 3. 

Based on the single-particle spectra of the ground state, one can construct excited states 
with high i^-values. In a deformed, axially symmetric nucleus, a high-i^ state is made by 
summing the contributions from several unpaired quasiparticles. To form low-lying high-fT 
states, several high-f2 single-particle (both neutron and proton) levels lying close to the Fermi 
surface are necessary. The well-deformed nuclei with A « 180, including 178 Hf, satisfy this 
requirement very well. With the restriction of the total Q and parity as 16 + , the candidate 
configurations of 178 Hf m2 will be constructed. 

For the ground state of 178 Hf, the neutron (proton) single-particle levels close to the 
Fermi surface are shown in the first column of the left (right) panel in Fig. [2j Each level is 
labeled by the Nilsson notation VL v [Nn z mi\ of its main component. It can be seen that in 
Fig. [21 the energy required to excite one neutron or one proton is not less than 0.63 or 1.84 
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MeV. As the experimental excitation energy of ^Hf" -2 equals 2.446 MeV, it is sufficient to 
consider one- or two-neutron excitations, together with one- or two-proton excitations. For 
two-neutron (two-proton) excitations, the following cases are considered: in the first case, a 
pair of particles below the Fermi surface are excited to two different levels above the Fermi 
surface; in the other case, two particles occupying different levels below the Fermi surface 
are excited to form a new pair. Those cases which involve four or more single-particle levels 
are not included for simplicity. All possible configurations are constructed by restricting the 
K (total Vt) value to 16 and the nuclear parity to +. Thus the configurations with the lowest 
excitation energies are obtained and labeled by 16 + ,16 + , 16 + , etc. 

The detailed configurations of the five lowest = 16 + states of 178 Hf are listed in column 
2 of Table [TJ In this table, the first 4 states are one- neutron plus one-proton excitation, while 
the last one is a pair of neutrons excited to high-f2 levels combined with the one-proton exci- 
tation. The configuration of 16+ is z/(7/2-[514])- 1 (9/2+[624]) 1 tt(7/2+[404])- 1 (9/2-[514]) 1 
with respect to the ground state, i.e., one formerly paired neutron in the level 7/2~[514] 
becomes unpaired and excited to the level 9/2 + [624], and another formerly paired proton 
in level 7/2 + [404] is excited to the level 9/2~[514]. This configuration is consistent with 
the former assignment of 178 Hf m2 l 2 ^^. Note that the excitation energy of 16 + given by 
the sum of single particle excitations is 3.954 MeV which is about 1.5 MeV higher than the 
experimental value 2.446 MeV. 

In order to obtain the self-consistent excitation energy for state 16 + microscopically, the 
diabatic constrained RMF calculations are carried out with the respective configuration 
information. As discussed before, the time-odd fields E 29 ^ caused by the unpaired nucleons 
should be taken into account carefully for the high-.fr state. The time-even calculation is 
also done for comparison. The PES of state 16 + with (without) current is plotted as a solid 
(dashed) red curve in Fig. [TJ together with the PES of state 16 + . The local minima of 16 + 
(16 + ) PES are denoted as up (down) triangles in Fig. [TJ For state 16 + , the excitation energy 
according to the calculations without and with current is, respectively 3.579 MeV and 2.801 
MeV, which clearly shows that both the self-consistent calculation and the consideration 
of the time-odd fields are crucial effects to obtain the reasonable excitation energy. The 
excitation energy of 16 + (E x = 2.801 MeV) is close to the experimental excitation energy 
(E x = 2.446 MeV) of ^Hf" 12 . The deformation (fii = 0.30) obtained for this isomer is similar 
to that of the ground state. In Table [TJ all the excitation energies from calculations without 
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and with current as well as the deformation calculated for the states 16J" to I65 are also 
given. For all five states, the excitation energies calculated with current are 0.6 ~0.8MeV 
smaller than those without current. By the way, the same approaches are also applied to the 
isomeric state 178 Hf m (experimentally K n =8~ ,E x =1.147MeV). The configuration obtained 
for 8^ is 7r(7/2 + [404]) _1 (9/2~[514]) 1 , and the excitation energy without and with current is, 
respectively 1.552 MeV and 1.315 MeV. 

In Fig. [2], the neutron and proton single-particle levels of state 16f from the RMF without 
and with current have been plotted in columns 2, 3 in both panels. In the time-odd calcula- 
tion each single-particle level splits into two levels due to the breaking of the time-reversal 
symmetry, the level with positive Q being energetically favored. Such splittings will change 
straightforward the energy gap between two single-particle levels. In particular, for one- 
neutron and one-proton excitation of the state 16* discussed here, the neutron energy gap 
between i/(9/2 + [624]) and z/(7/2 _ [514]) in the calculation without current decreases from 
2.00 MeV to 1.25 MeV , which is the gap between z/(9/2+[624], +Q) and i/(7/2 _ [514], -fi) 
in the calculation with current. Correspondingly the energy gap for the proton excitation 
concerned decreases from 1.34 MeV to 0.54 MeV. As a result, the calculation with current 
decreases the excitation energy considerably. 

In summary, the properties of the ground state of 178 Hf and the isomeric state ^Hf™ 2 
are investigated by the adiabatic and diabatic constrained RMF approaches. The con- 
strained RMF theory reproduces well the binding energy and deformation for the ground 
state of 178 Hf. With the single-particle levels of the ground state obtained in the adiabatic 
constrained RMF theory, by restricting K w = 16 + , the configuration with the lowest exci- 
tation energy is found to be z/(7/2-[514])- 1 (9/2+[624]) 1 tt(7/2+[404])- 1 (9/2-[514]) 1 , which 
is consistent with the former configuration assignment. The excitation energy based on the 
single-particle spectra is 3.954 MeV, which is much higher than the experimental excitation 
energy 2.446 MeV of 178 Hf m2 . By applying the self-consistent time-even and time-odd RMF 
calculation, the excitation energy of this configuration is decreased to 3.579 MeV and 2.801 
MeV, respectively. Therefore both the self-consistency and the consideration of current are 
important factors for studying nuclear isomers. 
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■ (adiabatic/diabatic) g.s. PES 
•16* PES (no /with current) 

•16* PES (no /with current) 




PK1 - 
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FIG. 1: (color online) Potential energy surfaces (PES) of Hf from adiabatic (open circles) and 
diabatic (solid curves) RMF calculations, as well as the PES of the first two 16 + states 16^ (red) 
and 16^" (blue) from the adiabatic constrained RMF calculations with (solid curves) and without 
(dashed curves) current. The RMF parameter set PK1 is adopted. The ground state as well as 
the local minima 16^~ (16^~) are denoted as a black asterisk and red (blue) triangles, respectively. 
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FIG. 2: (color online) Neutron (left panel) and proton (right panel) single-particle levels of the 
ground states (1st column), state 16^ of 178 Hf obtained from the constrained RMF approach 
without (2nd column) and with (3rd column) current. Each level is labeled by the Nilsson notation 
flir[Nn z mi] of its main component. The Fermi surfaces are marked by neutron and proton numbers 
in the grey circles. The red (blue) lines denote the parity + (— ), and the dashed (solid) lines denote 
the positive (negative) signs of O, respectively. The configuration difference is illustrated by filled 
circles, open circles and arrows. 
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TABLE I: The configurations, excitation energies as well as the deformations obtained from the 
adiabatic and diabatic constrained RMF approaches for the first five K n = 16+ states of 178 Hf. 
The excitation energies (in MeV) include the sum of single-particle excitations — and 

the excitation energies E x calculated without (with) current. 



State 


Configuration 




no 


E x 

current) (with 


current) (with 


02 

current) 


16 + 


tt(7/2+[404])- 1 (9/2-[514]) 1 


2.119 
1.835 


3.954 


3.579 


2.801 


0.302 


16+ 


i/(7/2+[633])- 1 (9/2+[624]) 1 
tt(7/2-[523])- 1 (9/2-[514]) 1 


1.549 
2.597 


4.146 


4.025 


3.342 


0.272 


16+ 


i/(5/2-[512])- 1 (11/2+[615]) 1 
7r(7/2+[404])^ 1 (9/2-[514]) 1 


3.204 
1.835 


5.039 


4.749 


4.129 


0.291 


16+ 


i/(7/2+[633])- 1 (9/2-[505]) 1 
tt(7/2+[404])- 1 (9/2-[514]) 1 


3.473 
1.835 


5.308 


5.380 


4.601 


0.273 


16f (5/2 


" [512])" 2 (9/2+ [624]) 1 (7/2- [503]) 1 
vr(7/2+[404])- 1 (9/2-[514]) 1 


4.113 
1.835 


5.948 


5.879 


5.302 


0.290 



12 



